home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b114.asm < prev    next >
Assembly Source File  |  1992-04-28  |  4KB  |  68 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.14    N Cascaded Real Biquad IIR Filters  
  8. ;   w(n) = x(n) - a1 * w(n-1) - a2 * w(n-2) 
  9. ;   y(n) = w(n) + b1 * w(n-1) + b2 * w(n-2) 
  10. ; 
  11. ;       X Memory Organization            Y Memory Organization 
  12. ; 
  13. ;                                                      b1N     Coef. + 4N-1 
  14. ;                                                      b2N     
  15. ;                                                      a1N     
  16. ;                                                      a2N     
  17. ;           wN(n-1)   Data + 2N-1                      .      
  18. ;           wN(n-2)                                     .      
  19. ;             .                                        b11     
  20. ;             .                                        b21     
  21. ;           w1(n-1)                                    a11   
  22. ;     R1,R0 ?  w1(n-2)   Data                    R4 ?    a21     Coef. 
  23. ;                                                          
  24. ; 
  25. ;
  26. ;                      DSP56000 IMPLEMENTATION 
  27. ;                                                           Program    ICycles 
  28. ;                                                           Words 
  29. ;       move  #$ffffffff,m0                                   2          2 
  30. ;       move  m0,m4                                           1          1 
  31. ;       move  #data,r0                                        2          2 
  32. ;       move  #coef,r4                                        2          2 
  33. ;       movep              x:input,a                          1          2 
  34. ;       move               x:(r0)+,x0    y:(r4)+,y0           1          1 
  35. ;       do     #n,end                                         2          3 
  36. ;       mac    -x0,y0,a    x:(r0)-,x1    y:(r4)+,y0           1          1 
  37. ;       macr   -x1,y0,a    x1,x:(r0)+    y:(r4)+,y0           1          1 
  38. ;       mac    x0,y0,a     a,x:(r0)+     y:(r4)+,y0           1          1 
  39. ;       mac    x1,y0,a     x:(r0)+,x0    y:(r4)+,y0           1          1 
  40. ;end 
  41. ;       rnd    a                                              1          1 
  42. ;       movep              a,x:output                         1          2 
  43. ;                                                 ------------------- 
  44. ;                                                      Totals 17       4N+16 
  45. ;                      DSP96002 IMPLEMENTATION 
  46.                                                              Program    ICycles
  47.                                                                     Words 
  48.        move   #$ffffffff,m0                                  ; 2     2 
  49.        move   m0,m4                                          ; 1     1 
  50.        move   m0,m1                                          ; 1     1 
  51.        move   #data,r0                                       ; 2     2 
  52.        move   r0,r1                                          ; 1     1 
  53.        move   #coef,r4                                       ; 2     2 
  54.        movep                       x:input,d0.s              ; 1     2 
  55.        fclr   d1                   x:(r0)+,d4.s y:(r4)+,d6.s ; 1     1 
  56.        do     #n,end                                         ; 2     3 
  57.        fmpy d4,d6,d1 fadd.s d1,d0  x:(r0)+,d5.s y:(r4)+,d6.s ; 1     1 
  58.        fmpy d5,d6,d1 fsub.s d1,d0  d5.s,x:(r1)+ y:(r4)+,d6.s ; 1     1 
  59.        fmpy d4,d6,d1 fsub.s d1,d0  x:(r0)+,d4.s y:(r4)+,d6.s ; 1     1 
  60.        fmpy d5,d6,d1 fadd.s d1,d0  d0.s,x:(r1)+ y:(r4)+,d6.s ; 1     1 
  61. end; 
  62.                      fadd.s d1,d0                            ; 1     1 
  63.        movep                       d0.s,x:output             ; 1     2 
  64. ;                                                            ; ---   --- 
  65. ;                                                     Totals:     19   4N+18 
  66. ;                                                            (17   4N+16) 
  67.  
  68.